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Abstract 

We investigate hadroproduction of B c mesons from initial gluons using QCD perturbation theory 
and nonrelativistic bound state approximations. Our results obtained by two completely indepen- 
dent calculations are confronted with existing results which contradict each other. Moreover, we 
examine the approximation based on heavy quark fragmentation and determine the range of validity 
of it. Predictions for cross sections and differential distributions at Tevatron and LHC energies are 
presented. 
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1 Introduction 



With the discovery of the top quark, it has become rather certain that B c mesons are the 
only flavored heavy quark resonances which exist in nature. Although they have not yet been 
observed experimentally, the production rates at the Tevatron may be sizeable enough to 
allow for experimental identification of these interesting bound states. Clear evidence can be 
expected from the LHC Jl|, where it should also be possible to measure specific production 
and decay properties. This will provide valuable tests of important aspects of strong and weak 
interactions. 

Because of the promising prospects, B c - physics is attracting increasing attention. In par- 
ticular, the hadronic production of B c mesons has been studied by several groups @-|| in 
the framework of QCD perturbation theory treating the B c couplings and wave functions in a 
nonrelativistic approximation. There is general agreement that in pp and pp collisions at high 
energies the gluon-gluon scattering process gg — > B c bc including higher order corrections is by 
far the dominant source of B c mesons. However, the quantitative predictions derived first in 
refs. H 0, |J in lowest order contradict each other in various respects. It is not always easy to 
compare the results, because of the use of different, sometimes outdated gluon momentum dis- 
tributions and the lack of clear statements on the choice of scales and bound state parameters. 
Yet, one can convince oneself that differences in the numerical input alone cannot explain the 
discrepancies. The problem is most obvious in refs. @ and Q which disagree already on the 
level of the subprocess. The confusion has also not been resolved in ref. [[| which appeared 
very recently while we were writing up this paper. The situation is rather unsatisfactory and 
needs clarification. 

Another controversial issue is the description of B c hadroproduction in terms of fe-pair produc- 
tion in gg-fusion followed by fragmentation of the 6-quark: gg — > bb, b — > B c bc. The applica- 
tion of this approximation [|[ was motivated by the success of the analogous approach to B c 
production in e + e~ annihilation. In fact, from the perturbative calculation of e + e~ — > B c bc 
[0| one can directly extract the relevant fragmentation function Di(z), z being the momen- 
tum fraction of the 6-quark transferred to the B c meson || ||. In hadronic processes, the 
hard scattering formalism involving convolutions of parton cross sections with structure and 
fragmentation functions is generally expected to work for high-p^ production. On the other 
hand, it has been shown in refs. JlO], [TTJ] that the fragmentation picture is not adequate for 
B c production in photon-photon scattering, not even at large 77 energies and the highest B c 
transverse momenta. Since the two processes 77 — > B c bc and gg — > B c bc are very similar at 
least if one disregards the gluon self-interactions and the color factors, one may mistrust the 
general expectation expressed above. While the previous investigations of this issue in refs. 
and dJ are not conclusive because of the discrepancies of these calculations, ref. [0 which 



also appeared just recently claims the failure of the fragmentation approximation. We will 
clarify this important question and show that the fragmentation description indeed works to 
the expected accuracy. 

In order to achieve these goals, we have performed two independent calculations of the dom- 
inant subprocess gg — > B c bc in 0(af) and in the usual nonrelativistic approximation. The 
result is compared numerically with the previous calculations. Where a comparison at the 
level of the subprocess was not possible, we have checked the convoluted cross sections for 
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the same parton distributions and input parameters as used in the calculation under exami- 
nation. For the integrated cross sections, we get partial agreement with ref. J|, however for 
a different value of mj, and only after including a factor 1/3. Although this missing factor 



has been corrected in ref. [12], there are still discrepancies, in particular, in p<r~distributions. 



Furthermore, we clearly disagree with refs. ||, |3|. || . 

For predictions on B c production at the Tevatron and LHC the gg cross section is folded with 



the updated gluon distribution of ref. [|T3]] in which the HERA measurements are taken into 
account. Moreover, we study in detail the dependence on the scale chosen for a s and the 
structure functions. Finally, we present distributions in transverse momentum and rapidity, 
as well as integrated cross sections including p T - and y-cuts. 

The paper is organized as follows. In chapter 2 we give details of the calculation of the 
matrix element and the phase space integration. There, we also specify the nominal values 
of all relevant parameters used in our numerical studies. Chapter 3 is devoted to total cross 
sections, and chapter 4 to differential distributions. A short summary and conclusions are 
given in chapter 5. 



2 Calculation of the subprocess gg —> B^b 



c 



To order a*, the subprocess gg — > B)*'bc is described by 36 Feynman diagrams. For the heavy 
h and c quarks, it is reasonable to neglect the relative momentum of the quark constituents 
and their binding energy relative to their masses. In this nonrelativistic limit, the constituents 
are on mass shell and move together with equal velocity. This implies the following relations 
for the masses and momenta of the c-quark, 6-quark, and the B^ bound state: 

M(B c )=M(B:) = M = m c + m b , Pc = j^P, Pi = -j±p. (1) 

Furthermore, the amplitudes for the production of S-wave states are obtained from the hard 
scattering amplitudes for gg — > bbcc by applying the projection operator JTJL ITT 



v(p- b )u(Pc) = -^0/ - M)Ii SSz , (2) 

where v{pi) and u(p c ) denote the respective quark spinors and Hss z is a s Pi n projector. For 
a pseudoscalar state one has n o = 75, while for a vector, Ilis z = j.. From now on, we 
differentiate between the pseudoscalar ground state B c and the vector ground state B*. The 
decay constants appearing in eq. (0) are defined by the matrix elements 

< 0\h^c\B c (p) > = if Bc p„, < Q\h,c\B:(p) > = M/ fl .e M , (3) 
and are related to the bound state wave function at the origin by 

/B o = /flS = ^l*(0)|. (4) 

As in ref. the color structure is not accounted for in eq. (0). The color factors are 

computed separately when squaring the amplitude. 
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In order to have an internal cross-check, we have calculated the matrix element squared by 
two independent methods. One calculation is based on the usual trace technique, while the 
other calculation makes use of helicity amplitudes similarly as in ref. [lTf. The latter method 



is a generalization of the technique described in detail in ref. [15] for bosonic final states in 
e + e~ scattering. The results of these two calculations are compared numerically for different 
sets of particle momenta and found to agree to seven digits. Furthermore, it has been checked 
that the matrix element vanishes when the polarization vector of anyone of the initial gluons 
is substituted by the momentum vector of this gluon. 

The phase space integration is also performed independently using two different routines. 
In a first step, the integration routines are checked for the pure 3-particle phase space and 
found to agree numerically with the corresponding analytical formulae to four digits for all 
distributions considered. Before integrating the matrix element squared with the help of the 
Monte Carlo routine VEGAS fpT6|| , the strongest peaks of it were smoothed out by introducing 
new integration variables. After the integration all differential distributions and cross sections 
of the two independent calculations agree within the Monte Carlo errors. 

In addition, we have calculated cross sections and distributions for bb(cc) production and 
subsequent fragmentation of the 6(c)-quark, that is for gg — > bb, b — > B c c and gg — > cc, c — > 
B c b, respectively. The relevant fragmentation functions D b (z) and D c (z) are known from 
perturbation theory ||. In e + e~ — > B c cb, they provide a perfect approximation for the energy 
distribution da/dz with an error of the order of M 2 /s, y/s being the cm. energy [|. In 
contrast, in 77 — ► B c cb above threshold, the fragmentation mechanisms are dominated by 
quark recombination processes. This can be shown in a gauge invariant way without actually 



factorizing the process in heavy quark production and fragmentation |pT| . Because of the 
similarity of the underlying processes, one anticipates recombination dominance in gg — ► B c cb 
at least at low p?- Concerning the fragmentation mechanism, we stress that the factorized 
description da(gg — > bb) £g> D b (z) is not justified at small px and close to threshold, where the 



quark masses are important. This caveat is ignored in ref. [JR]. In order to find out where 
and to what extent this approximation can be trusted, we have also assumed various relations 
between daughter and parent momenta: 



p T = s/A - p? sin 9 b , p L = p T cot 6 b} \i = M (I), m b (II), (III), (5) 

where is the gluon-gluon cm. energy. The case I respects the physical phase space 
boundaries, II and III are assumed in refs. || and ||12|| , respectively. These prescriptions are 



equivalent to the accuracy of the fragmentation approach. 

To obtain cross sections for pp or pp collisions, the gluon-gluon cross section a has to be folded 
with the gluon structure functions g(x,Q 2 ) of the (anti)proton: 



o is) 



[ dxi [ dx 2 g(xi,Q 2 )g(x2,Q 2 )9(xix 2 s - AM 2 )a(sx 1 x 2 ,Q 2 )- (6) 

JO JO 



Analogous convolution formulae hold for differential distributions. For numerical illustrations 
we use the following input 

f Bc = f B = 0.4 GeV, m b = 4.8 GeV, m c = 1.5 GeV (7) 
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together with the MRS (A') gluon distribution [[13]. The strong coupling constant a s {Q 2 ) 
entering a, is taken in leading logarithmic approximation: 



a Ami 



-i . 33— 2nr 

1 + -E»r a s 



ln(- 



(8) 



with rtj = 5 and a s (m 2 z ) = 0.113. This is the normalization of the running coupling in the 
MRS(A ; ) parton distributions. As is usually done in such calculations, we identify the scale of 
a s with the evolution scale of the structure functions. 




200 400 600 800 1000 

s 1/2 , GeV 

Fig. 1. Integrated cross sections for gg — > 
versus the cm. energy. The lower (up- 
per) curve corresponds to B c (B*) production. 




0.001 1 1 — 1 — 1 — ' 

10 100 1000 10000 

s 1/2 , GeV 

Fig. 2. Total cross sections for pp(pp) — > 
(*) 

B y c 'bc + X versus the cm. energy. The lower 
(upper) curve corresponds to B c (B*) produc- 
tion calculated with the scale Q 2 = p\ + M 2 . 
The symbols refer to other choices. 



3 Total cross sections 



First we present the parton cross sections a(s) for the subprocess gg — > By 'be, and the 
convoluted production cross sections for pp( and pp) — > B^bc + X. Then we confront our 
results with results in the published and unpublished literature 0, |3], [|, |5| . 

Figures 1 and 2 illustrate our main findings. The partonic cross sections are shown in Fig. 1 for 
a fixed value of a s in order to facilitate comparison with other calculations. The production 
threshold is located at \fl = 2(m& + m c ) = 12.6 GeV. At sufficiently high energies, the cross 
sections fall roughly like 1/s. Moreover, the B* cross section is approximately three times as 
large as the B c cross section as expected from spin counting. 

The total hadronic cross sections for B c and B* production are plotted in Fig. 2. One can 
see the typical rise of a with energy reflecting the rise of the gluon density as x approaches 
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20 GeV 


40 GeV 


100 GeV 


1 TeV 




8.27(2) 


12.7(1) 


6.67(4) 


0.208(2) 




29.3(1) 


40(1) 


21(1) 


0.27(4) 




20.3(1) 


33.1(2) 


17.4(1) 


0.61(1) 




71.6(3) 


105(2) 


56(2) 


0.6(1) 



Table 1: Comparison of total cross sections in pb for gg — > B^bc with the corresponding 
results of ref. [4]. The input parameters are nib = 5.1 GeV, m c = 1.5 GeV, a s = 0.2 and 
Jb c = Jb* = 0.57 GeV. The number in parenthesis shows the Monte Carlo uncertainty in the 
last digit. 



%min — 4(m& + m c ) 2 /s, and the peaking of a near threshold. A similar behaviour is observed 
for production in the scattering of two bremsstrahlung photons radiated from e ± -beams 
[ [TI|| . The scale dependence is also exhibited in Fig. 2 for a variety of choices including the 
fixed scale Q 2 = AM 2 . At relatively low energies, the predictions vary by more than one order 
of magnitude. At higher energies, however, the results become more stable changing only 
within a factor two. It is interesting to note in this respect that the differences resulting from 
convolution with the gluon densities MRS(A ; ) of ref. [Tj| and CTEQ2 of ref. [|T7| are very small. 



In fact, they would be invisible in Fig. 2 . In addition to the scale uncertainty which can only 
be reduced by calculating the next-to-leading corrections, a formidable task indeed, one has 
uncertainties connected with the decay constants f B (*) and the effective quark masses nib and 
m c . These amount at least to another factor of two. 

Next, we confront our results with the numerical results of refs. || ||, |], The comparison is 
easiest and clearest for the parton cross section a, where one does not have to pay attention to 
different parametrizations of the gluon distribution and for a constant value of a s . Needless to 
say, we have always adopted the same values for nib, m c and / as the reference calculation. 

Table [I] summarizes the comparison with ref. ||. As can be seen, we clearly disagree with 
this calculation. Amazingly, if one uses instead of the quark masses quoted the masses given 
in eq. (|?|) one achieves agreement at = 20, 40 and 100 GeV within the Monte Carlo errors, 
up to a color factor 1/3 which has apparently been missed. In the recent repetition l!| of 



their calculation, the authors account for the missing color factor. However, we do not see an 
explanation for the remaining discrepancy at V§ = 1 TeV which exceeds the MC uncertainty 
by far. We also observe a clear discrepancy with ref. || as demonstrated in Table 2. The fact 
that these authors claim rough agreement with ref. speaks against their results because of 
the error in the color factor pointed out above. 

Furthermore, results on the gg-cross section can be found in the second paper of ref. |2j. 
Unfortunately, since the choice of the QCD-scale A and of the momentum scale Q 2 is left 
unclear, a direct comparison is difficult. In Table 0, we try two hypotheses: fixed a s and 
a s (es), where the numerical factor e can be chosen arbitrarily. As one can see, in none of the 
two cases we are able to achieve agreement. 

Finally, the hadronic cross sections given in ref. || seem to differ even by more than one order 
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20 GeV 


30 GeV 


40 GeV 


60 GeV 


80 GeV 


100 GeV 




6.86(2) 


9.71(4) 


9.74(5) 


7.93(5) 


6.23(5) 


5.01(4) 




23.6 


26.0 


22.3 


15.5 


11.2 


8.4 



Table 2: Comparison of the total cross section in pb for gg — > B c bc with the result of ref. [6]. 
The input parameters are m h + m c = 6.3 GeV, m& = 3m c , a s = 0.2 and f Bc = 0.50 GeV. The 
number in parenthesis shows the Monte Carlo uncertainty of the last digit. 



Vs [GeV] 


Ref. @ 


fixed a s 


a s (es) 


20 


4.9 


4.9 


4.9 


30 


8.5 


7.0 


5.0 


60 


7.9 


5.9 


2.4 



Table 3: Comparison of the total cross section a(gg — > B c bc) (in pb) for fixed and running a s 
with the result quoted in ref. [2]. The input parameters are = 4.9 GeV, m c = 1.5 GeV, 
and fs c = 0.48 GeV. The normalization of a s is chosen such as to reproduce the cross section 
at 20 GeV. 



of magnitude from our results. Again, since the normalization of a s (Q 2 ) is not specified, we 
cannot make a more definite statement. 

To conclude, all previous calculations |2], |3], ffl, Q are in mutual disagreement. Our calculation 
confirms the results on a of ref. [[| at < 100 GeV when including the missing color 
factor ]12| and taking m;, = 4.8 GeV. However, at ^/l = 1 TeV, a clear numerical discrepancy 
remains. 



4 Distributions 



Having discussed the total production cross sections, we now turn to the differential 
distributions in transverse momentum and rapidity. Since they are very similar for B c and 
B*, we show them for the pseudoscalar meson only. Moreover, since the sensitivity to the 
scale in a s (Q 2 ) and g(x,Q 2 ) is already indicated in Fig. 2, we shall present here results only 
for the scale Q 2 = p\ + M 2 . Finally, we will examine the fragmentation description and 
determine the region of validity of this approximation. In order to be consistent with the 
complete calculation at the fixed order a*, we do not take into account evolution effects in 
the fragmentation function Di(z). These effects are studied in ref. [0. Moreover, it should be 
mentioned that c-quark production and fragmentation is negligible. 

Fig. 3 illustrates the main features of the subprocess gg — > B c bc. We see that the fragmentation 
description approaches the ^^distributions resulting from the complete calculation only in 
the tails of the distributions. At \f§ < 40 GeV a problem arises also when px approaches p™ ax 
because of mass ambiguities in the phase space boundaries of the fragmentation approach. 
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These effects increase as s decreases and are largest for the prescription III of eq. (|5|). On the 
other hand, if the physical phase space boundary is imposed on the fragmentation approach 
as in eq. (|5]) case I, the approximation slightly improves. Comparing Fig. 3 with Fig. 5 of ref. 
II], one sees that in contrast to 77 — > B c bc where heavy quark fragmentation is completely 
subdominant even at large pt, in gg — > B c bc it provides a good approximation at large s and 
Pt- Obviously, the presence of gluon self-couplings and color factors has a drastic influence on 
the relative importance of the fragmentation and recombination. 

dcr/dp T , pb/GeV 



0.001 





. • 

00, 


a s = 0.2 ; 


- c\ • 
o\ H 




o\ • 




: °o\ 

cm 




\^ u\ 


I — 


° 

/III ods^ Q 



dcr/dp T , pb/GeV 
1000 



20 30 40 /s \ , 50 

p T , GeV 




Fig. 3. Transverse momentum distributions of 
the B c in gg — ► B c bc at VI = 40 and 100 GeV: 
complete 0(ct^) calculation (circles) and frag- 
mentation approximation (solid curves). The 
labels I-III refer to the kinematics specified in 
eq. (5). 



Fig. 4. Transverse momentum distribution of 
the B c at Tevatron and LHC energies: com- 
plete 0(a^) calculation (circles) and fragmen- 
tation approximation (solid curves). The labels 
I and III refer to the kinematics specified in eq. 
(5). At each energy, results are shown without 
a rapidity cut (full circles) and for \y\ < 0.5 
(empty circles). 



Distributions in px for pp (and pp) collisions are displayed in Fig. 4. Also shown is the 
effect of a stringent cut in rapidity. A few comments are in order. Most importantly, after 
convolution the full calculation and the fragmentation description are in reasonable agreement 
at pt > 10 GeV. This can be understood from the properties of the unfolded j>r~distributions 
illustrated in Fig. 3 and from the rise of the gluon density at small x which favours contributions 
from the smallest possible subenergies s. At pr < 10 GeV, the deviation grows indicating 
that one is getting outside the range of validity of the fragmentation picture. Furthermore, 
the sensitivity to the kinematical prescription, eq. (|5|), decreases slowly with increasing pt- 
Beyond the p<r~ r egion considered in Fig. 4, it can be safely ignored. These observations are 
rather independent of the rapidity cut. The rapidity distributions exhibit the usual plateau 
region around y = 0, as shown in Fig. 5. 



Thus, in contradiction to ref. fl~2 |, we find that the fragmentation description provides a 
good approximation for p T > 10 GeV. The deviation from the complete 0(af) results is much 
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smaller than the uncertainties in the scale Q 2 , in the decay constant /„<*), and in the effective 
values of mj and m c . 



dcr/dy, pb 

10000 h 



100 



0.01 



Q 2 = p T 2 + M 2 



-10 



=1 4TeV 




/ min\ 

a(p T >p T ), pb 

10000 



100 



1 r 



0.01 




Fig. 5. The rapidity distributions of the B c 
mesons at Tevatron and LHC energies. 



Fig. 6. Integrated production cross sections 
for the B c as functions of the minimum px 
cut: complete 0(a 4 ) calculation (circles) and 
fragmentation approximation using eq. (5) II 
(solid curves). For each energy, results are 
shown without a rapidity cut (full circles) and 
for \y\ < 1 (empty circles). 



5 Summary and conclusions 

We have presented new results on B c meson production in pp and pp collisions derived from 
the 0(a 4 ) subprocess gg — > B^'bc using nonrelativistic bound state approximations. We 
have performed two independent calculations employing two different methods. We therefore 
believe that our results which deviate from all previous calculations @, [|. [|. |12|] are correct. 

In addition, we have determined the range of validity of the fragmentation approximation 
with the heavy quark fragmentation function Di(z) as calculated in perturbation theory M. 
For pt > 10 GeV, we find good quantitative agreement between the fragmentation description 
and the complete calculation. At pt = 5 — 10 GeV, mass ambiguities in the fragmentation 
kinematics lead to effects up to a factor of two. Below 5 GeV, the fragmentation approximation 
brakes down. 

In order to evaluate the observability of B c mesons at the Tevatron and LHC, it is useful 
to integrate the ^^distributions of Fig. 4 over p? > pjp m . The resulting cross sections are 
plotted in Fig. 6. This estimate indicates a production of the order of 10 4 B c mesons with 
pmtn _ Q e y a £ tn e Tevatron, assuming an integrated luminosity of 100 pb -1 . Contributions 
from the production and decay of B* mesons and heavier states are not yet taken into account 
in this number (see e.g. ref. ||). This rate may be sufficient for a first observation in the 
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decay channels B c — > J/tpX, for which a branching ratio of the order of 10% is predicted |18| . 
Finally, at the LHC for lOOffcT 1 one can expect 10 7 direct B c mesons at p™ n = 20 GeV which 
should make B c physics a very interesting topic. 
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